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We propose a model for the photo-deactivation mechanism for DNA based upon accurate quantum 
chemical and molecular dynamical evaluations of model Watson/ Crick nucleoside pairs and stacked 
pairs. Our results corroborate recent ultrafast experimental studies on DNA oligonucleotides and 
suggest that following photo-excitation to a local it — n* state, the excitation is rapidly delocalized 
over several (3-4) bases on an ultrafast time-scale. However, this delocalized state is unstable with 
respect to the motions of the protons involved in hydrogen-bonding between Watson/Crick pairs 
and rapidly re-localizes to a charge-transfer state on a longer time-scale ranging from 10 to 100 
ps. This state, too, is unstable and relaxes via a conical intersection with the ground state near 
the geometry of the enol- and imino-tautomeric form. We suggest that this internal deactivation 
mechanism is responsible for the intrinsic photostability of DNA. 



The mutagenic effect of the UV irradiation on the ge- 
netic code has been a topic of interest in the last few 
decades. As a result of the absorption of the short- 
wavelength component of UV light (A < 290 nm) by 
the stratospheric ozone layer the spectrum of sunlight 
at the sea level barely overlap the low energy absorp- 
tion tail of the nucleic acid bases: adenine (A) , thymine 
(T), guanine (G), and cytosine (C). However, the de- 
pletion of the ozone layer in the recent years raises the 
concern of increasing exposure to the UV radiation. Ex- 
isting reports indicate that for every 1% decrease in the 
ozone layer there could be 4% increase in the incidence 
of skin cancer. Photoexciation of DNA in the UVB 
and UVA region gives rise to photolcsions, primarily in 
the form of cyclobutanc pyrimidine dimcrs (CPD) and 
the pyrimidine(6-4)pyrimidone photoadducts. [TJ [2] The 
presence of these lesions are associated with a number 
of nonmelanomic skin cancers. [3] The adverse effect of 
these photolesions is largely mitigated via enzymatic ac- 
tion involving DNA photolyase; however, such secondary 
mechanisms are both error prone and energetically costly 
and most certainly developed later in the evolutionary 
process. Experimental evidence also suggests that pho- 
toexcitations in DNA can be rapidly quenched through 
internal conversion and dissipated as heat. [4, 5 

Until development of the femtosecond laser spec- 
troscopy, attempts to study the dynamics of these ex- 
cited states were thwarted by the ultrafast character of 
the internal conversion process. In the last decade the 
combined effort of the researchers employing femtosecond 
spectroscopy, molecular beam techniques, and quantum 
chemistry methods shed light on the excited-state dy- 
namics of the nucleic acids. However the precursor state 
which initiates the photochemical reactions damaging the 
genome remains elusive. Ultrafast spectroscopic stud- 
ies by Fiebig et al. [51 H] indicate that following photo- 
excitation, the exciton rapidly delocalizes over a number 
of stacked bases. A similar explanation of the long-lived 
excited states was offered by Markovitsi and coworkers. 
They interpret the results of the fluorescence upconver- 
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FIG. 1: Ab initio ground and excited state potential curves 
for a dAdT Watson/Crick nucleoside pair along the N6(A)- 
H stretching coordinate. At each point along the curve, the 
ground-state geometry was optimized constraining the sugars 
to their positions in a B-DNA chain. 

sion and time-correlated single photon counting experi- 
ments in the frame of the exciton theory. In this model 
the large number of neutral excitons states are formed 
over several bases located both on the same strand and 
on opposite strands with subsequent relaxation of the ex- 
citon. In this approach both base stacking and base pair- 
ing determine the excited state dynamics. [8 ., 9 . ITUl [TT1 IT^] 
Similarly, by comparing the ultrafast transient absorp- 
tion signals for isolated bases, single-stranded oligonucle- 
osides and base-paired oligonucleosides, Kohler's group 
concludes that because the ground state recovery of sin- 
gle and double strand systems are essentially identical 
and very different from isolated bases. [T3j [14j [15j [TBI [17] 
These transient absorption on dAdT oligonucleotides 
support a model whereby the initial local excitation 
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evolves on an ultrafast time scale to form a long-lived 
species that subsequntly decays to the electronic ground 
state. Kohler et al. assigns the long-lived species to 
inter-base charge-transfer states in which an electron is 
transferred from one base to a 7r-stacked neighbor to form 
an intrastrand charge-transfer exciplex. While the pre- 
cise nature of these long-lived states has been subject to 
debate, [17j HE] it is clear from these results that base 
stacking interactions as opposed to the electron-driven 
proton-transfer dictates the excited state dynamics in 
double-stranded DNA. 

In contrast to these results, experiments and calcula- 
tions by Domcke and coworkers [HI [2Qj |2T] on isolated 
base-pairs and pyridine dimers indicate that upon pho- 
toexcitation to a irir* state, the transfer of a proton from 
one molecule to the other carries the system through two 
conical intersections, first between an excitonic tttt* state 
and wtt* charge-transfer (CT) state and second between 
the CT state and the ground-state on a timescale of sw 
60-100 ps in matrix-isolated model base-pair dimer sys- 
tems depending upon system. Based upon these results, 
a back and forth proton transfer could account for the 
intrinsic photo-stability in actual DNA chains. 

Here, we report on our series of theoretical investiga- 
tions aimed at providing a unified dynamical model for 
describing the internal photochemical relaxation path- 
ways in a DNA oligomer. [22] [23] Our studies combine 
accurate and high-level quantum chemical descriptions 
of the relevant excited states of a Watson-Crick AT nu- 
cleoside pair. Stacking interactions are introduced by 
computing the exact Coulomb coupling between excited 
state transition densities on neighboring base-pairs as 
sampled from a molecular dynamics simulation of a 12- 
base DNA oligomer of d(A)-d(T) in water in its "B" form. 
By including both base-stacking, base-pairing, and pro- 
ton exchange on an equal footing, we hope to provide 
an unbiased assessment of the initial photophysical pro- 
cesses. Our results indicate that photoexcited states in 
DNA chain can be rapidly localized by fluctuations in the 
protons involved with base pairing and that excited-state 
proton transfer between paired nucleosides can effectively 
quench photoexcitations. Weak coupling between de- 
localized neutral excitonic states and localized charge- 
transfer states does account for the enhanced lifetimes of 
the photoexciations. 

The chemical structure of the dAdT base-pair is shown 
in the inset of Fig. [l] along with several key bond dis- 
tances. In agreement with previous calculations both 
equilibrium structures were found slightly buckled and 
propeller- twisted. [3D] In addition, the distance between 
the two fragments, A and T, is shorter in the "rare" tau- 
tomer, a feature also calculated by Villani.[24] Particu- 
larly, the N6-H-04 hydrogen bridge in the imino/enol- 
tautomer is shorter by 0.324 A compared with the cor- 
responding distance in the Watson/Crick structure. The 
energy profile shown in Fig. [T|was constructed by mov- 
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FIG. 2: Plot of the oscillator strength corresponding to the 
nir* and nn* states as a function N6(A)-H bond distance. 
The difference densities shown next to each trace correspond 
to the ground state configuration of the WC AT base pair. 



ing the N6 proton of adenine toward thymine and reop- 
timizing all other bases geometrical parameters. These 
states are energetically very close together and resolving 
them accurately requires both large electronic bases-sets, 
accurate treatment of electronic-correlations, and the in- 
clusion of virtual excitations. Furthermore the influence 
of these effects is magnified as the system is pushed to- 
wards the two conical intersections encountered as we 
stretch the rN6—H bond. Finally, we found that there is 
is considerable interplay between the bases and the sug- 
ars that we could not simply replace them with methyl 
groups. The barrier for the double proton transfer dAdT 
— > dA*dT* in the ground state calculated at the HF/cc- 
pVDZ is rather high, 31.1 kcal/mol, while the barrier for 
the reverse process dA*dT* — > dAdT calculated at the 
same level amounts to 14.9 kcal/mol (Fig|lJ. Especially, 
the latter barrier is significantly higher compared with 
the values of 0.2 and 0.02 kcal/mol reported by Gorb 
et al[25] and Guallar et al. |26j respectively, for a AT 
base pair without sugars. However, reoptimization of 
the structures with sugars at the MP2/def-SV(P) level 
lowers the barriers for the forward and reverse processes 
to 18.2 and 2.0 kcal/mol respectively (data not shown). 

Vertical excitation energies (VEE) of the lowest mr* , 
localy excited 7T7t* excitonic states and charge transfer 
7T7T* states of the Watson-Crick base pair calculated as a 
function of the N6-H distance of adenine at the CIS(D) 
level are plotted in Fig. [I] At the ground state equi- 
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FIG. 3: The difference densities of stacked 9H-adenine pen- 
tamer (left) calculated at the CIS(D)/cc-pVDZ level corre- 
sponding to two different localied nn* states (middle and 
right. 



librium geometry, the calculated lowest energy excited 
state corresponds to the mr* state with both the n and 
ir* orbitals completely localized on thymine. The 7T7t* 
exciton states are just about 0.2-0.4 eV above the mr* 
state. We also determine that state corresponding to the 
charge-transfer (tttt* CT) state is approximately 0.6 eV 
above the ittt* localized exciton states. The vertical exci- 
tation energy calculated at CIS(D) level are generally in 
good agreement with those calculated at the CC2 level 
for the AT base pair without sugars using similar quality 
basis set [27]. The primary configurations in all three 7T7t* 
exciton states correspond to electron/hole excitations lo- 
calized on cither thymine or adenine bases. The weights 
of these configurations are by and large similar, below 
50%, producing the irir* excited states delocalized over 
both bases with none of them entirely localized on just 
one base. The calculated vertical excitation energy of the 
7T7T* exciton states increase while the vertical excitation 
of the 7T7r* CT state decreases with increasing N6-H (ade- 
nine) bond distance up to 1.4 A. As the adenine N6-H 
bond distance becomes longer this trend is reversed. Fur- 
thermore, the character of the excited states also changes. 
The 7T7r* exciton states become increasingly localized on 
just one base cither adenine or thymine with the con- 
tribution from the corresponding excitation greater than 
80%. Ultimately, this is related to the formation of the 
tautomer, When the dAdT base pair geometry is reopti- 
mized with the N6-H (adenine) bond distance r# = 1.5 
A the N3 hydrogen of thymine moves over to the other 
side forming a covalent bond with the adenine Nl to form 
the imino/enol tautomer. 

We also computed the excitation energies for stacked 
pcntamer of 9H-adcnine and saw little evidence of charge 
transfer exciplex states between stacked bases. None of 
the 20 excited states calculated at the CIS(D)/cc-pVDZ 
level have a charge transfer character. The lowest energy 



states can be most adequately classified as mr* and irir* 
states. Likewise the vertical excitation energies (VEE) of 
the AT base pair the VEEs of the former states of stacked 
pcntamer arc slightly lower than the VEEs of the latter 
states. The irir* states with the largest oscillator strength 
are also the most delocalized, the corresponding molecu- 
lar orbitals being delocalized over all five adenine bases. 
On the contrary, the mr* states can be localized on a 
single thymine base. Preliminary CIS(D) calculations on 
stacked adenosine pentamers based upon MD configura- 
tions give little evidence for exciplex formation, although 
depending upon the instantaneous geometry of the chain, 
excitons may have a small static dipole due to differing 
degrees of localization of the occupied and virtual orbitals 
contributing to the configuration-interaction expansion 
of the excited state wave functions. The difference densi- 
ties corresponding to the 7T7t* states of the pentamer with 
the largest oscillator strength are shown in Fig. [3j These 
calculations were performed in vacuum, consequently, it 
is entirely possible that contributions from solvent po- 
larization, counter-charges, and near by water molecules 
could stabilize intrastrand exciplex states. One can argue 
against exciplex formation in stacked homodimers since 

the exciton binding energy (sa 0.4 0.5 eV typical for 

conjugated organic species) is far greater than the differ- 
ence in either electron affinities or ionization potentials of 
the stacked pair. [2H] [23] Based upon the exciplex stabil- 
ity criteria, exciplex states are only expected for stacked 
heterodimers (eg. AT, CG, etc..) 

For the purpose of developing a minimal model for the 
excited states of DNA chain we need to make a number 
of judicious assumptions. The Coulombic coupling ele- 
ments between two stacked Watson-Crick AT base pairs 
were calculated using the transition density cube (TDC) 
method described in our previous report [23], The tran- 
sition densities have been calculated for several Watson- 
Crick AT base pairs molecular geometries extracted from 
the molecular dynamics simulations. In Fig. [2] we show 
the transition densities calculated at the CIS/cc-pVDZ 
level corresponding to two different localied 7T7t* states 
superimposed at the neighboring AT base pairs of the 
(dA) 10 .(dT) 10 . 

Our quantum chemical results suggest that a minimal 
model for the relevant states of a DNA oligomer must 
consist of a ground state \g) and localized mr* excitons 
( | Ej ) ) and local 7T7t* charge-transfer states (\cj)) for a 
each Watson-Crick pair. We assume that only these 
states are involved in the derealization and relaxation 
process following photoexcitation. We assume that the 
overlap between neighboring states is sufficiently small 
that (e n \e n >) = (c„|cv) = 5 nn > and (e„|c„) = 0. The cou- 
pling interactions between neighboring excitonic states 
|e,-) and |ej±i) (A) were determined to be in the range 
200-400 cm -1 . This sets the time scale for coherent ex- 
citon transfer between stacked bases to be on the order 
of 5-10 fs. 
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FIG. 4: The three-dimensional potential energy surface de- 
scribing the motion of protons between N6(A) and 04(T) 
and between N3(T) and N1(A) shows two critical points 
in the ground state. The deeper minimum corresponds to 
the amine/keto structure of AT and a shallow one to the 
imine/enol structure (A*T*). Upon absorption of a UV pho- 
ton (blue arrow) the initaly delocalized excitonic states (1) 
undergo a rapid localization on ~10 ps timescale for single 
bases and «100 ps timescale for stacked base pairs to form 
a charge transfer (CT) states (2). The subsequent CT states 
passing through a conical intersection are carried back to the 
ground state (3). 



Fig. [4] shows the adiabatic potentials arising from our 
reduced model. First, the lowest surface is the poten- 
tial surface for tautomerization in the ground state. Two 
minima occur at (ri,r2) = (3,0) and (0,3) correspond- 
ing to one base pair or the other in the tautomer form. 
However, both excitonic states are unstable with respect 
to the proton exchange coordinates. Once the system 
has moved away from the origin along one of the proton- 
transfer coordinates, the electronic states rapidly localize 
and we are carried towards the conical intersection be- 
tween the local CT state and the ground state. 

Let us assume that the lifetime of the delocalized state 
is limited by proton transfer between one of the base 
pairs such that as soon as one proton coordinate cross 



The excitonic couplings are generally much stronger 
than those between the local excitonic and CT states (j'2 
= 50 cm' 1 ) or between two stacked CT irir* states (\cj) 
and |cj±i)) (20 cm -1 ) as calculated using the transition 
density cube. These stacking interaction are highly sen- 
sitive to the relative orientation of base pairs and can 
change by an order of magnitude. However, the absolute 
magnitude of the interaction is very small compared with 
the vertical excitation energies of the d(A)-d(T) base pair 
calculated at CIS(D)/cc-pVDZ level. Consequently, if we 
estimate the energy lowering of excitonic tttt* states due 
to the excitonic coupling as twice the nearest-neighbor 
coupling calculated using TDC method the delocalized 
exciton is about 800 cm -1 lower in energy than a local- 
ized exciton. Finally we include coupling between the 
the ground state \g) and the 7T7t* states \ej) (ji) and the 
ground state and the CT states (j'3). We systematically 
adjusted these terms as to best reproduce the energies of 
the ground state and the lowest energy 7T7t* states of the 
d(A)-d(T) nucleoside from our ab initio studies. (FigjlJ. 
The adiabatic energy surfaces along the proton-transter 
normal mode coordinate r n cx (6(a)N6-H + $(t)N3-h) 
for a dimer can be obtained from the eigenvalues of the 
following 5x5 matrix: 



(1) 



I 

the XT/CT intersection, the delocalized state collapses 
to from a localized state. Assuming the usual Condon 
separation between the nuclear and electronic dynamics, 
we can write this within the non-adiabatic Marcus ap- 
proximation 

kloc = ^ |2 1 =e -(AE +Er )»/(4E r fc B T) (2) 

where Vif = j% is the diabatic coupling. We can esti- 
mate this rate by setting the driving force AE to be the 
energy difference between the vertical exciton and the 
ground-state tautomer and the reorganization energy E r 
as E r = E ct — E t . This sets the time-scale for interbase 
electron transfer of r = l/fc/ oc = 10.3 ps. This gives a 
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lower limit to the exciton lifetime since even a small er- 
ror in our values can change this by a factor of 2 to 5. 
Moreover, for the delocalized case, the coupling matrix 
element will be at least proportional to the probability for 
finding the exciton on a given site, V a b oc jiPn- Thus, for 
the delocalized case where the exciton is extended over 
3-5 bases we expect r w 100 to 250 ps. In summary, 
this study combined with our previous study of exciton 
derealization in B-DNA chains [23 proposes the follow- 
ing mechanism. Following vertical n — n* excitation of 
an adenosine, the exciton rapidly delocalizes between 3- 
4 neighboring stacked A's on a time-scale given by the 
exciton-exciton coupling, A. The derealization length is 
limited by the fact that A is strongly modulated by the 
structural fluctuations of the DNA chain about its ideal 
B-DNA form. This initial derealization occurs on the 
femtosecond time-scale. Next, these states are unstable 
with respect to the fluctuations of the stretching motions 
of the protons involved in the hydrogen-bonding between 
the Watson/Crick pairs. This causes a re-localization of 
the exciton to occur on the time-scale of pico-seconds. 
Subsequent relaxation to the ground-state occurs on a 
longer time scale as determined by the conical intersec- 
tion between the CT and the ground-state. 

Although this study has focused solely on poly(dA)- 
poly(dT) chains, the implication for more general se- 
quences is clear. It is entirely possible that such in- 
ternal photodeactivation mechanisms are responsible for 
the intrinsic photo-stability of DNA. Such mechanisms 
may have provided a necessary evolutionary advantage 
for early DNA-based life forms. 
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